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Abstract 

VT (= parity times time-reversal) symmetry of complex Hamiltonians with real 
spectra is usually interpreted as a weaker mathematical substitute for Hermiticity. 
Perhaps an equally important role is played by the related strengthened analyticity 
assumptions. In a constructive illustration we complexify a few potentials solvable 
only in s— wave. Then we continue their domain from semi-axis to the whole axis 
and get the new exactly solvable models. Their energies come out real as expected. 
The new one- dimensional spectra themselves differ quite significantly from their 
s— wave predecessors. 
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1 Introduction 



Our mathematical understanding of many physical systems can become drastically 
simplified after their suitable complexification. This is true, first of all, in the study 
of resonances [1] and of several other quantum scattering phenomena [2]. Recently, 
the idea of working in a complexified phase space for bound states [3] re-entered the 
scene with a new enthusiasm supported by an immediate relevance of the related 
break-down of parity V in certain field theories [4]. 

In the mathematically more accessible quantum mechanical models certain 
exceptional complex interactions with VT symmetry happen to become strictly 
equivalent to a real potential after a supersymmetric [5] or integral, Fourier-like [6] 
transformation. For other models, the analysis of the related purely real spectra 
of energies has been performed by several techniques. One may recollect, e.g., the 
most straightforward numerical experiments [7] , semiclassical approximants [8] and 
the so called delta expansions [9]. Resummations of divergent perturbation series 
[10] and the so called exact WKB method [11] also offered several Hamiltonians 
for which the spectra of energies E n were proved strictly real. 

One of the most immediate sources of information about the possible con- 
nection or correlation between the absence of a decay Im E n = and the VT 
symmetry H = VTHVT itself is provided by the exactly solvable models in one 
dimension. Step by step, there were proposed the VT symmetric versions of the 
harmonic oscillator [8], of the asymmetric Morse interaction [12], of the asymp- 
totically symmetric (sometimes called "scarf ) hyperbolic oscillator [13] and of its 
asymptotically asymmetric but locally not too dissimilar (also known as Rosen- 
Morse) alternative [14]. Before complexification, all of them belong among the 
so called shape invariant potentials (cf. the review [15]) so that some of their 
properties can be clarified using the language of supersymmetry [16]. 

On the basis of numerical experience [9] the current attention is exclusively paid 
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to the forces V(x) which are analytic in x. An extremely interesting byproduct of 
this point of view can be found in a transition to more dimensions for quartic (i.e., 
unsolvable) oscillators [6] and for the central and exactly solvable VT symmetrized 
harmonic oscillator [17] and Coulomb problem [18]. Within the set of the similar 
forces with a centrifugal-like singularity there still exist a few models without a clear 
interpretation. After a glimpse in the Table 4.1 of the review [15] we immediately 
discover two of them, namely, the Eckart model 

V^\r) = A{A - l) -2B C ^ (1) 
sinh r sinhr 

and the generalized Poschl- Teller potential 

cosh r sinh r 

In the standard interpretation [2], both these s— wave models are only partially, 
incompletely solvable and, in this sense, lie somewhere in a "territory of nobody" . 
This was the main source of our present inspiration. We see no reason why these 
two interactions should not be appropriately continued to the whole line and clas- 
sified, afterwards, as the two new or "forgotten" exactly solvable VT symmetric 
models. In detail, this will be done in sections 2 and 3 below. 

We have to remind the reader that the later force (2) may be often found in the 
current literature in its alternative form V^ GPT \r) — (u + v cosh 2r)/sinh 2 2r [15] 
or in the special form known as Hulthen potential [2] . The former correspondence 
is mediated by the trivial re-scaling of the axis of coordinates by factor 2. In the 
present context, the latter, much less trivial relationship deserves a more explicit 
attention. Its thorough discussion will be added here, therefore, in section 4. A 
few further relevant overall comments may be found in our summary and final 
discussion in section 5. 
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2 VT— regularization and Eckart oscillator 



From the purely historical point of view the loss of Hermiticity in the domain 
of complex couplings proved more than compensated by the new insight in the 
solutions of one of the most popular unsolvable models V(x) = uox 2 + Ax 4 [19]. 
Today, its spectrum is understood as a single multi-sheeted analytic function of the 
complex coupling constant A 6(T. The same idea applies to the set of resonances in 
the cubic well V(x) = uj x 2 + A x 3 although a careful analytic continuation must be 
also performed in the coordinate x itself [7]. These observations guided the semi- 
classical and numerical studies of the forces V^(x) = u x 2 + g x 2 {i x) 5 containing a 
variable real exponent 5 [9]. The related VT— symmetric quantum mechanics with 
its new perturbation series [20] as well as quasi-classical approximation schemes 
[11] and matrix-truncation methods [21] works with the globally, asymptotically 
deformed paths of integration in the related Schrodinger equation 

if>(x) = Eif>(x). (3) 

The 5 = 2, quartic anharmonic oscillator of ref. [22] exemplifies these systems 
which need not remain integrable on the real line. Its asymptotic integrability and 
decrease of wave functions is only recovered after we bend both our coordinate 
semi-axes downwards and replace 

{x>l} — > {x = Qe- iv }, {x<-l} — > {x = -Qe iv } 

beyond certain distance g 1 and within certain bounds upon ip e (0, 7r/3). The 
further growth of S beyond 5 = 2 would make both the asymptotical if— wedges 
shrink and rotate downwards in the complex plane. 

Let us defer the discussion of the similar cases to our last section 4 below. Re- 
turning now just to our first two examples (1) and (2) we may notice that both of 
them may be characterized by a "weak" , ip = option. Globally they do not leave 



3 



the real axis of x at all. Such a simplification proves most natural in the 5 — > oo 
regular model of ref. [23]), admitting the most natural physical interpretation of 
the real physical coordinates after all. The related VT— symmetrized oscillators 
need not necessarily differ from their Hermitian counterparts too much. One can 
hope to encounter just slight modifications of the formulae available, e.g., in the 
factorization context [24] and in its Lie-algebraic [25], operator [26] or supersym- 
metric [27] re-interpretations. 

Equally straightforward innovations may be expected in the domain of our 
singular forces (1) and (2). One can simply avoid their isolated singularities by a 
local deformation of the integration path. In this way the strong repulsion in the 
origin (so popular in some phenomenological models [28] and fully impenetrable in 
one dimension) becomes readily tractable via a suitable choice of the cut. 

2.1 Terminating solutions revisited 

Once we pay attention to the real s— wave potential (1) with the strongly singular 
core, usually attributed to Eckart [29], we have to keep in mind that this Hermitian 
model is solvable on the half-line only, with r G (0, oo) and, conventionally, A > 1/2 
and B > A 2 . Its fixed value of the angular momentum i = is in effect a non- 
locality which lowers its practical relevance in three and more dimensions. 

As already mentioned, the local deformation of the integration path enables us 
to forget about the strong singularity in the origin. We may admit the presence of 
the so called irregular components in ip(r) ~ r 1 - A near r = 0. They would be, of 
course, unphysical in the usual formalism [30]. Here, on the contrary, we continue 
r — > x with x G (— oo, oo) and encounter the new possibilities. 

In the new perspective we have to re-analyze the whole Schrodinger equation 
anew. Our choice of appropriate variables 

= to- !)•(, + DV(^), y = saf = i-* 
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is dictated by the arguments of Levai [27], and the consequent VT — symmetry 
considerations require that we use the purely imaginary couplings B = if3. Then 
we insert V^ Eck \x) in eq. (3) and our change of variables leads to its new form 

z(l - z) <p"(z) + [c - (a + b + l)z] ip'(z) - abip(z) = (4) 

where 

c=l + 2w, a + b = 2u + 2v + l, ab = (u + v)(u + v + 1) + A(l - A) (5) 
and 

Av 2 = 2B-E, Au 2 = -2B - E. (6) 

Our differential equation is of the Gauss hypergeometric type and its general so- 
lution is well known [31], 

<p(z) = C x ■ 2 F 1 {a, b; c; z)+C 2 - z 1 ^ 2 F 1 (a + 1 - c, b + 1 - c; 2 - c; z). (7) 

Besides the obvious relevance of such an exceptional solvability of a model with 
a strong singularity in quantum mechanics, an independent encouragement of its 
study is provided by its methodical appeal in the context of field theory, especially 
in connection with the so called Klauder phenomenon [32]. 

2.2 Asymptotic boundary conditions 

Technically, the first thing we notice is that our parameters a and b are merely 
functions of the sum u + v and vice versa, u + v = (a + b — l)/2. The immediate 
insertion then gives the rule (a - b) 2 = (2A - l) 2 and we may eliminate 

a = b±(2A-l). (8) 

We assume that our solutions obey the standard oscillation theorems [33] and 
become compatible with the boundary conditions ^(±oo) = in eq. (3) at a 



5 



discrete set of energies, i.e., if and only if the infinite series 2-^1 terminate. Due 
to the complete a <-> b symmetry, we only have to distinguish between the two 
possible choices of C 2 = and C± — 0. 

In the former case with the convenient b = —N (= non-positive integer) the 
resulting numbers a + b and u + v prove both real. Using the definition of B the 
difference u—v = —i/3/(u+v) comes out purely imaginary. The related terminating 
wave function series (7), i.e., 

/ 1 \ u+v 

*<*> = (i£hi) ^ < 9 > 

is asymptotically normalizable if and only if u + v > 0. This condition fixes the 
sign in eq. (8) and gives the explicit values of all the necessary parameters, 

a = 2A - N -1, u + v = A-N-l, u-v = -i- A — ^ . (10) 

A — N — I v ' 

For all the non-negative integers N < N max < A — 1 the spectrum of energies is 
obtained in the following closed form, 

E = -\ (u 2 + v 2 ) = -(A-N-l) 2 + {A _^_ ir N = 0,l,...,N max . 

(11) 

The normalizable wave functions become proportional to Jacobi polynomials, 

<p[z(x)] = const. ■ Pjv /2 ' v/2) (cothx). (12) 

We have shortly to return to the second option with C\ = in eq. (7). Curiously 
enough, this does not bring us anything new. Although the second Gauss series 
terminates at the different b = c — 1 — N, the factor z l ~ c changes the asymptotics 
and one only reproduces the former solution. All the differences prove purely 
formal. In the language of our formulae one just replaces u by — u in (and only in) 
both equations (9) and (10). No change occurs in the polynomial (12). 
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3 Poschl- Teller potential 

Schrodinger equation (3) with the bell-shaped potential V(r) ~ l/cosh 2 r belongs 
to the most popular exactly solvable models in quantum mechanics. Its appli- 
cations range from the analyses of stability and quantization of solitons [34] to 
phenomenological studies in atomic and molecular physics [35], chemistry [36], 
biophysics [37] and astrophysics [38]. Its appeal involves the solvability by differ- 
ent methods [27] as well as a remarkable role in the scattering [2]. Its bound-state 
wave functions represented by Jacobi polynomials are also encountered as super- 
partners of a complex "scarf model [13]. 

Not too surprisingly, virtually all these applications lose their physical ground 
after an addition of the repulsive spike. Still, it is not too difficult to extend the 
exact solvability itself to the latter potential called, often, the Poschl- Teller well 
[39]. The related Schrodinger equation (3) must be confined to semi-axis r e (0, oo) 
or appropriately regularized. 

3.1 Regular izat ion 

We may repeat that the impossibility of using the real V^ PT ^ of eq. (2) with 
A > B > in more dimensions nor on the whole axis in one dimension is felt 
unfortunate in methodical considerations and in perturbation theory [40]. Singu- 
larities of the centrifugal type are encountered in phenomenological models [28, 41] 
but, unfortunately, not too many of them are solvable [42]. 

In our present regularization of the singularity we shall not deform the straight 
integration path at all. We shall rather proceed in a way inspired by the pioneering 
paper [6] where Buslaev and Grecchi employed simply a constant downward shift 
of the whole coordinate axis, 

r = x — ie, x e (—00,00). (13) 
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In a way similar to the oscillator V^ BB \x) = V^ HO \x—ic) = x 2 —2icx—c 2 of ref. [8] 
and to its three-dimensional generalization [17] the meaning of the VT symmetry 
degenerates here to the mere trivial invariance with respect to the simultaneous 
reflection x — > —x and complex conjugation i — * —i. The shift (13) is the main 
source of a regularization here. As long as l/(x — is) 2 — (x + ie) 2 /(x 2 + e 2 ) 2 at 
any e ^ 0, the centrifugal term remains nicely bounded in a way which is uniform 
with respect to x. Without any difficulties one is able to work with the similar 
centrifugal-like terms on the whole real line of x. 

The same idea applies to the regularized Poschl- Teller potential 

V^ RPT \x) = V^ PT \x -is), 0<e< tt/2. 

This potential is a simple function of the Levai's [27] variable g(r) = cosh2r. As 
long as g(x — ie) = cosh 2x cos2£ — i sinh2x sin 25, the new force is VT symmetric 
on the real line of x G (— oo, oo), 

V { - RPT \-x) = [V { - RPT \x)Y. 

Due to the estimates | sinh 2 (a;— is) \ 2 = sinh 2 x cos 2 e+cosh 2 x sin 2 e = sinh 2 x+sin 2 e 
and | cosh 2 (a; — ie)\ 2 = sinh 2 a; + cos 2 e the regularity of V^ RPT \x) is guaranteed 
for any parameter e G (0, tt/2). 

3.2 Solutions 

In a way paralleling the the preceding section the mere analytic continuation of the 
s— wave bound states does not give the complete solution. One must return to the 
original differential equation (3). There we may conveniently fixA+l/2 = a>0 
and B -1/2 = (3 > and write 

( d 2 (3 2 -l/A a 2 -l/A\ , . . ^ , . . 

~ + 7- 2 ' -yf-Ui =^^. r(x)=x-te. (14) 

\ dx z smh r(x) cosh r(x) ) 
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This is the Gauss differential equation 



z(l + z) <p"(z) + [c + (a + b + l)z] tp'(z) + a& <p(z) = 



(15) 



in the new variables 



(1 + ^(4 



z 




using the suitable re-parameterizations 



o? = {2v-\j2)\ 



P 2 = (2/i 



1/2) 2 , 



2/i + l/2 = c, 2fi + 2u = a + b, 



E = -{a-bf. 



In the new notation we have the wave functions 



i){x) = smh^ +1 iy{x)}cos\f a+1 ' 2 [r{x)} V [z{x)} 



(16) 



with the sign ambiguities r = ±1 and a — ±1 in 2/i — t/3+1/2 and 2v = aa+ 1/2. 
This formula contains the general solution of hypergeometric eq. (15), 

<p(z) = Ci 2 F 1 {a, b- c; -z) + 2 Fi(a + 1 - c,b + 1 - c;2 - c; -z). (17) 

The solution obeys the complex version of the Sturm-Liouville oscillation theorem 
[33]. In the case of the discrete spectrum this means that we have to demand 
the termination of our infinite hypergeometric series, suppressing its undesirable 
asymptotic growth at x — > ±oo. 

In a deeper analysis let us first put C 2 — 0. We may satisfy the termination 
condition by the non-positive integer choice of b = —N. This implies that a = 
N + 1 + era + t/3 is real and that our wave function may be made asymptotically 
(exponentially) vanishing under certain conditions. Inspection of the formula (16) 
recovers that the boundary condition ^(±oo) = will be satisfied if and only if 



1 < 2N + 1 < 2N max + 1 < -aa - t(3. 
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The closed Jacobi polynomial representation of the wave functions follows easily, 

The final insertions of parameters define the spectrum of energies, 

E = -(27V + 1 +aa + r{3) 2 < 0. (18) 

Now we have to return to eq. (17) once more. A careful analysis of the other 
possibility C\ — does not recover anything new. The same solution is obtained, 
with r replaced by — r. We may keep C2 = and mark the two independent 
solutions by the sign r. Once we define the maximal integers which are 

compatible with the inequality 

2N%£ + l<-aa-r(3 (19) 

we get the constraint iV < N^J . The set of our main quantum numbers is finite. 

4 Bent contours and Hulthen potentials 

In both our above examples (1) and (2) an overall VT symmetry of the Hamiltonian 
is, presumably, responsible for the existence of the real and discrete spectrum [8]. 
Cannata et al [16] and Bender et al [43] were probably the first to notice that one of 
the various limits 5 — > 00 of the power-law models with ip — > n/2—0(l/S) becomes, 
unexpectedly, exactly solvable again, in terms of special Bessel functions. These 
observations attract attention to strongly deformed contours. One of possibilities 
of their interpretation is the Liouvillean change of variables [44] . 
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4.1 The VT symmetry preserving changes of variables 

In the first step let us recollect that in the spirit of the old Liouville's paper [45] 
the change of the (real) coordinates (say, r <-> £) in Schrodinger equation 

d 2 



dr 



2 + W(r) 



x (r) = -k x(r) 



(20) 



mediates a transition to a different potential. In terms of an invertible function 



r = r(£) which possesses a few first derivatives 


r'(0, r 




. we g< 


bound state problem with the new interaction 








V(0-E=lr>(0] 2 {WHO] + K 2 } + l 


\r"(0~ 
[r'(0_ 


2 1 

2 


\r"'(0] 
r'(£) 



(21) 



and normalizable wave functions 



(22) 



In the Jacobi-polynomial context the Liouvillean changes of variables have been 
applied systematically to all the Hermitian models (cf. Figure 5.1 in the review 
[15] or ref. [46] for a more detailed illustration). A similar exhaustive study is still 
missing for the VT symmetric models within the same class. Let us now try to 
partially fill the gap. For the sake of brevity we shall only restrict our attention to 
the VT symmetric initial eq. (20) with the Poschl- Teller potential 

(3 2 - 1/4 a 2 - 1/4 



W(r) = 



x — is, 



x G (— oo, oo). (23) 



sinh 2 r cosh 2 r ' 
Its normalizable bound states are proportional to the Jacobi polynomials, 

X (r) = sinh r/3+1/2 r cosh CTQ+1/2 r P^ aa) (cosh 2r) 

at all the negative energies — k 2 < such that 

K = k { ^ t) = -aa - t/3 - 2n - 1 > 0. 
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These bound states are numbered by n = 0, 1, . . . ,n^J and by the generalized 
parities a — ±1 and r = ±1. 

We may note that our initial VT symmetric model (20) remains manifestly 
regular provided only that its constant downward shift of the coordinates r = 
r \x) = x — ie remains constrained to a finite interval, e G (0, 7r/2). In a key step of 
its present modification let us now change the coordinates as follows, 

sinhr = —ie 1 ^, £ = v — iu. (24) 

This shifts and removes the singularity at r = to infinity (u — > +oo). In an 
opposite direction, one cannot proceed equally easily from a choice of a realistic 
V(£) to the re-constructed coordinate r(£). This methodical asymmetry is due to 
the definition (21) containing the third derivatives and, hence, too complicated to 
solve. Still we are quite lucky with our purely trial and error choice of eq. (24). 
Firstly, the real line of x becomes mapped upon a manifestly VT symmetric curve 
£ = v — iu in accordance with the compact and invertible trigonometric rules 

sinhxcose = e u sinv, 
cosh x sine = e u cosv, 

i.e., in such a way that 

v = arctan (^f ) = v {x) E («(_«,), «(«,)) = (-f + e, f - e) , 
u = U( x ) = | In (sinh 2 x + sin 2 e^j . 

Our path of £ is a down-bent arch which starts in its left imaginary minus infinity, 
ends in its right imaginary minus infinity while its top lies at x — v — and 
— u = — M(o) = In 1/ sine > 0. The top may move towards the singularity in a 
way mimicked by the diminishing shift e — > 0. Although the singularity originally 
occurred at the finite value r — > 0, it has now been removed upwards, i.e., in the 
direction of — u — > +oo. 
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4.2 Consequences 

The first consequence of our particular change of variables (24) is that it does not 
change the asymptotics of the wave functions. As long as r'(£) = ?tanhr(£) the 
transition from eq. (20) to (3) introduces just an inessential phase factor in 
This implies that the normalizability (at a physical energy) as well as its violations 
(off the discrete spectrum) are both in a one-to-one correspondence. 

The explicit relation between the old and new energies and couplings is not too 
complicated. Patient computations reveal its closed form. With a bit of luck, the 
solution proves non-numerical. The new form of the potential and of its binding 
energies is derived by the mere insertion in eq. (21), 

"'«) = (r^ + r^. E = « 2 - < 25 > 

At the imaginary £ and vanishing A = this interaction coincides with the Hulthen 
potential. 

In the new formula for the energies one has to notice their positivity. This is 
extremely interesting since the potential itself is asymptotically vanishing at both 
ends of its integration path. One may immediately recollect that a similar paradox 
has already been observed in a few other VT symmetric models with an asymptotic 
decrease of the potential to minus infinity [22, 47]. 

The exact solvability of our modified Hulthen potential is not yet guaranteed 
at all. A critical point is that the new couplings depend on the old energies and, 
hence, on the discrete quantum numbers n, a and r in principle. This could induce 
an undesirable state-dependence into our new potential. Vice versa, the closed 
solvability of the constraint which forbids this state-dependence will be equivalent 
to the solvability at last. A removal of this obstacle means in effect a transfer of 
the state-dependence (i.e., of the n— , a— and r— dependence) in 

A = A(a) = 1 - a 2 , C (= A + B) = k 2 - f3 2 
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from C to (3. To this end, employing the known explicit form of k we may re-write 



C = C(a, r, n) = (era + 2n + l)(cra + 2n + 1 + 2t/3). 



(26) 



This formula is linear in t/3 and, hence, its inversion is easy and defines the desirable 
state-dependent quantity (3 = f3(a,T,n) as an elementary function of the constant 
C. The new energy spectrum acquires the closed form 



E = E(a,T,n) = A + B + 



1 



era + 2n + 1 — 



A + B 



aa + 2n + 1 



(27) 



Our construction is complete. The range of the quantum numbers n, a and r 
remains the same as above. 



5 Discussion 



5.1 Spectrum of the VT symmetric Eckart model 

The new spectrum of energies seems phenomenologically appealing. The separate 
N— th energy remains negative if and only if the imaginary coupling stays suffi- 
ciently weak, f3 2 < (A — N — l) 4 . Vice versa, the highest energies may become 
positive, with E = E(N max ) growing extremely quickly whenever the value of the 
coupling A approaches its integer lower estimate 1 + N max from above. In this 
way, even a weak VT symmetric force V^ Eck \x) is able to produce a high- lying 
normalizable excitation. This feature does not seem connected to the presence of 
the singularity as it closely parallels the similar phenomenon observed for the VT 
symmetric Rosen- Morse oscillator which remains regular in the origin [14]. Also, 
in a way resembling harmonic oscillators the distance of levels in our model is 
safely bounded from below. Abbreviating D = A — N — 1 = A e ff ective > its easy 
estimate 

£ "- £ ~-> = < 2D + 1 >( 1 + D^W) >1 
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(useful, say, in perturbative considerations) may readily be improved to — 
E N _i > (3 2 /D 2 at small D < 1, to E N - E N ^ > 2D at large D > 1 and, in 
general, to an algebraic precise estimate obtainable, say, via MAPLE [48]. 

Let us emphasize in the conclusion that the formulae we obtained are completely 
different from the usual Hermitian s— wave results as derived, say, by Levai [27]. 
He had to start from the regularity in the origin which implied an opposite sign in 
eq. (8). This had to end up with the constraint B > 0. Moreover, the size of B 
limited the number of bound states. 

In the present VT symmetric setting, a few paradoxes emerge in this com- 
parison. Some of them may be directly related to the repulsive real core in our 
V {Eck \x) with imaginary B. Thus, one may notice that the increase of the real 
repulsion lowers the iV— th energy. In connection with that, the number of levels 
grows with the increase of coupling A. In effect, the new bound-state levels emerge 
as decreasing from the positive infinity (!). At the same time, the presence of 
the imaginary B = if3 shifts the whole spectrum upwards precisely in the manner 
known from non-singular models. 

5.2 Paradoxes in the Poschl Teller case 

Let us now compare our final result (18) with the known e = formulae for s waves 
[27]. An additional physical boundary condition must be imposed in the latter 
singular limit. This condition fixes the unique pair a — — 1 and r = +1. Thus, the 
set of the s— wave energy levels En is not empty if and only if a— f3 > 1. In contrast, 
all our e > potentials acquire a uniform bound \V^ RPT {x)\ < const < oo. Due 
to their regularity, no additional constraint is needed. Our new spectrum E^' T ^ 
becomes richer. For the sufficiently strong couplings it proves composed of the 
three separate parts, 

4r- } <0, 0<N<N^~\ a + (3>l, 
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E ( N ' +) <0, 0<N<N!~£\ a>(3 + l, (28) 

< 0, < N < N£>-\ 0>a+l. 

The former one is non-empty at A + B > 1 (with our above separate conventions 
A > —1/2 and B > 1/2). Concerning the latter two alternative sets, they may 
exist either at A > B or at B > A + 2, respectively. We may summarize that 
in a parallel to the VT symmetrized harmonic oscillator of ref. [17] we have the 
Nmdx^ + 1 quasi-odd or "perturbed" , analytically continued s— wave states (with 
a nodal zero near the origin) complemented by certain additional solutions. 

In the first failure of a complete analogy the number N^~"> + 1 of our quasi- 
even states proves systematically higher than N^~^ + 1, especially at the larger 
"repulsion" > 1. This is a certain paradox, strengthened by the existence of 
another quasi-odd family which behaves very non-perturbatively. Its members 
(with the ground state ipQ + ' \x) = cosh A+1 [r(x)]sinh 1 ^ B [r(x)] etc) do not seem 
to have any s— wave analogue. They are formed at the prevalent repulsion B > 
A + 2 which is even more counter-intuitive. The exact solvability of our example 
enables us to understand this apparent paradox clearly. In a way characteristic for 
many VT symmetric systems some of the states are bound by an antisymmetric 
imaginary well. A successful description of its perturbative forms V(x) = ux 2 + 
i\x 3 [7, 10] carries numerous analogies with the real and symmetric V(x) = ux 2 + 
Ax 4 . The similar mechanism creates the states with (a, r) = (+, — ) in the present 
example. 

A significant novelty of our new model V^ RPT ^ (x) lies in the dominance of its 

imaginary component at the short distances, x ~ 0. Indeed, we may expand our 

force to the first order in the small e > 0. This gives the approximation 

1 = smh 2 (x + ie) = 1 + ^ cosh a; + Q ^ ^ 

sinh 2 (x — is) (sinh 2 x + sin 2 e) 2 sinh 2 x sinh 3 x 

We see immediately the clear prevalence of the imaginary part at the short dis- 
tances, especially at all the negligible A = 0(e 2 ). 
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An alternative approach to the above paradox may be mediated by a sudden 
transition from the domain of a small e ~ to the opposite extreme with e ~ tt/2. 
This is a shift which changes cosh a; into sinhx and vice versa. It intertwines the 
role of a and (3 as a strength of the smooth attraction and of the singular repulsion, 
respectively. The perturbative/non-perturbative interpretation of both our quasi- 
odd subsets of states becomes mutually interchanged near both the extremes of 
the parameter e. 

The dominant part (29) of our present model leaves its asymptotics compara- 
tively irrelevant. In contrast to many other VT symmetric models as available in 
the current literature our potential vanishes asymptotically, 

An introduction and analysis of continuous spectra in the VT symmetric quantum 
mechanics seems rendered possible at positive energies. This question will be left 
open here. 

In the same spirit we may also touch the problem of the possible breakdown 
of the VT symmetry. In our present solvable example the violation of the VT 
symmetry is easily mimicked by the complex choice of the couplings a and (3. 
Due to our closed formulae the energies will still stay real, provided only that 
Im (a a + t(3) = 0. 

5.3 Transition to the Hulthen model 

In the light of our new results we may now split the whole family of the exactly 
solvable VT symmetric models which contain a strong singularity in the two dis- 
tinct categories. The first one "lives" on the real line and may be represented or 
illustrated not only by the popular Laguerre-solvable harmonic oscillator [17] but 
also by both our present Jacobi-solvable forces. The second category requires a 
arch-shaped path of integration which lies confined within a narrow vertical strip. 
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It also involves both the Laguerre and Jacobi solvable subsets. The former one 
may be represented by the complex Morse model of ref. [12] and by the Coulomb 
force with a complex charge [18]. Our present new Hulthen example offers their 
first Jacobi-solvable counterpart. The parallels may be illustrated by the following 
picture 



symmetric 
Vl H °)(r) [17] 

X r = —i exp i x 



periodic 
V^ M \x) [12] 



symmetric 
V( PT )(r) [49] 



| sinhr = — iexpix 




where the vertical correspondence originates from the changes of variables. One 
notices the similarities in the (symmetric or periodic) form of the functions V as 
well as the differences in the straight-line or bent-curve shapes of the domains 
r — r(t) E(D or x = x(t) G(Z7, respectively. 

The less formal difference between the two categories may be also sought in 
their immediate physical relevance. Applications of the former class may be facil- 
itated by a limiting transition which is able to return them back on the usual real 
line. In contrast, the second category may rather find its most useful place in the 
methodical considerations concerning, e.g., field theories and the mechanisms of 
the parity breaking [4]. Within the quantum mechanics itself the second category 
might also parallel the studies of the "smoothed" square wells in non-Hermitian 
setting [16, 43]. 

In the conclusion let us recollect that the VT symmetry of a Hamiltonian 
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replaces and, in a way, generalizes its usual Hermiticity. This is the main reason 
why there exists an unexplored space for new solvable models. In their context, 
an example with an "intermediate" , hyperbola-shaped arc of coordinates remains 
still to be discovered. Up to now this type of contour has only been encountered 
in the "quasi-solvable" (i.e., partially numerical) model of ref. [22]. 
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